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NATIONAL ADVISORY COMMH?TEE FOE AERONAUTICS 


TECHNICAL MEMORANDUM NO. 1215 


FLOW PATTERN IN A CONVERGING-DIVERGING NOZZLE* 
By El. Oswatltsch and W. Rothstein 


The present report describes a new method for the prediction 
of the flow pattern of a. gas in the two-dimensional and axially 
symmetrical case. It Is assumed that the expansion of the gas is 
adiabatic and the flow stationary. The several assumptions 
necessary on the nozzle shape effect, in general, no essential 
limitation on the conventional nozzles. The method Is applicable 
throughout the entire speed range,* the velocity of Bound itself 
plays no singular part. The principal -reight is placed on the 
treatment of the flow near the throat of a converging-diverging 
nozzle. For slender nozzles formulas are derxvod for the calcula- 
tion of - the velocity components as function of the location. 


I. INTRODUCTION 


The field of a compressible nozzle flow has been treated 
repeatedly. Thus, Th. Meyer (reference l) computed the. transition 
from subsonic to supersonic 'flow for a given velocity distribution 
over the nozzle axis. G. I. Taylor (reference 2) calculated the 
case of subsonic flow wtdch at the throat of the nozzle reaches 
such speeds that the velocity of sound is exceeded at Beveral 
points. H. Gortler (reference 3) dealt in particular with the 
transition from one of these types of flow into the other. 

H. DERIVATION OF THE FUNDAMENTAL EQUATIONS 
FOR TWO-DIMENSIONAL FLOW 


The nozzle axis is indicated with x, the normal to it is y, 
and the origin of the coordinate system is .placed in the center of 
the narrowest cross section of a Laval nozzle symmetrical about 
the x-axis, figure 1. The velocity components are u and v, the 
direction of flow Is frcm left to right. The shape of tho nozzle 
Is given by a function f(x), with f denoting half the height 

*"Das Stromnngsf eld . in oiner Lavalduse , u Jahrbuch 1942 
Luftfahrtforschung, pp. I 91-102. 
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of the nozzle. Taking f(o) = 1 all linear dimensions ■with half 
the smallest nozzle height become nondimensional. 

In contrast to the earlier reports u- and v are expressed 
by a power series in y, tlie coefficients of which depend on x. 
In view of the nozzle shape being, symmetrical to the x-axis the 
flow itself is visualized as symmetrical to the x-axiB. The 
velocity at the x-axis is denoted by the subscript 0. Vq in 
this case is zero. Limited to terms of the fourth power of y 
the velocity components are ■ • • 


u(x, y) = Uq + agy 2 + Jy + . . . 
v(x, y) = t x y + ~y b 3 y 3 + . . . 


* 




(l) 


with Uq ; & 2 j a^j b^ and as fundtions of x. The contents of 

the present report consist in establishing the relationship 
between these quantities and function ■ f(x). The first task will 
be to ascertain the relationship of the coefficients eg, a],, 
b-,, and. b^ with the quantities Uq(x) and f(x). 


(a) Calculation of the Coefficients ag, . b^, and bg 


Since a potential flow is always required,- the following 
equation of irrotational motion 


gives by means of (l) 


du c3v 
oy c3x 


( 2 ) 



h 1 'y + 



The prime mark on b^ 1 and b^* is to -indicate the derivation with 

respect to x. Between the coefficients the' following relations- 
must therefore prevail: • * 


a s - h 



(2a) 
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The density p is made dimensionless by the “tank density, " that 

is, the density at speed v = yu^ + v^ = 0 and ell velocities "by 
the maximum velocity, that is, the velocity at density p = 0. 
Thus the energy theorem takes the form 


T? - 1 - p 14 - 1 (3) 

K is the ratio of the specific heats .» 

Elimination of pressure and density frcm the Euler equation 
of the continuity condition and the adiabatic equation, leaves the 
dynamic- gas equation 


( 0 2 . u2) |a + (o2 . h: . 2ot *».. 0 

oy ay / • 


ox 


cr is the squared velocity of sound, which may be -written as 


c 2 „ iLz_i pK-l . llJ: (i _ 

2 2 K ' 


ib) 


Elimination of the sonic velocity from the dynamic-gas equation 
by means of this equation gives 



cftr 

k - l / 5y 


(■ 


K + 1 

>r“~i 


u 





3u 

c)y 


(5) 


Forming by (l) the first derivatives of u and v and entering 
these along with the velocity components themselves in equation (5) 
gives by comparison of coefficients the relations between uo, 

a 2 > a lp hp hg, ', a^ and Uq'. The terms independent, 
of y give a particularly simple result. As is readily apparent 



1 - u 0 2 



(6) 


With equation (2a) the quantity a 2 in its relation with u 0 
and their derivatives can be computed. Furthermore, it is pointed 
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out at this point that the terms bjy and 


-l a 2? i 2 

2 


in (l) relative 


to u q are comparatively small for a nozzle flow. Thus for a 
first approximation of a 2 and b-j_ it is sufficient tosinsert the 

velocity afforded hy the simple flow filament theory and its 
derivative in equation (6). 


With u 0 as the flow filament velocity the, equation of 
continuity for the simple flow filament theory reads 


1 

u s (^l - Ug 2 ^" 1 f = const. 

A. 

In compressible flow the function u^l - u B 2 ) K_i substitutes for 
the velocity, for incompressible' flow the quantity u s . A special 
symbol for this function is found to be practical. We put 

1 

,e(u) i u(i - u 2 ) K_1 ... ( 7 ) 

The variation of the function is for K = 1,400" as follows, the 
values of the functions are given in tabl 9 I. 

For the derivatives of’ 9 with respect to u the following 
abbreviations are introduced 


i k + 1 _ g 

1 $3 _ 0 u ' k - T d _ l d% 6 uu 
"e du * e = T 1^2 = ~ 

and 

G ■ 

2 



(7a) 


(7b) 


If the function, for the velocity at' the axis u^ ■ of for the 

flow filament velocity. Ug is needed, it is indicated by the 
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subscript o or s. The formulas for b^ and ag can then be 
•written as 



The terms 


3 and — a^y 1 *' play a prominent part at the 


3! b 3 y “ — im 

nozzle edge, since they disappear -with the third end fourth power 
of y with approach to the x-exis. Hence it is logical to use 
these terms to satisfy the boundary conditions. With u^ and vg 


designating the Telocity components at the nozzle edge, the following 
boundary conditions must be fulfilled: 


= VLgl 


Ft 


(9) 


or also 


V + |V 3 = f '(»0 + l a 2 f2 + iii b 3' fl ‘) 

Hence a differential equation for b^ which can be written in the 
form 


| b 3 f3 = f ' (u 0 4- I a 2 f 2 ) - V 4. -A- (3c) 

The term not underscored is small compared to the term f *Uq, as 

seen from (l). It probably can be scored as a rule. At least it 
makes equation (8c) easy to iterate, by scoring the last term at 
the firBt step of iteration and then, by graphical differentiation 
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of the t >2 obtained in first approximation, calculating a second 

approximation of this quantity. With the nearly formed derivative 
this b^ ' scarcely differs from the previous one. It should be 

noted that it is not a matter of the relative variation of b n 

1 j, 

but that the variation of the term — b„ 'f 1- in relation to u„ 

2k 3 0 

is decisive. This iteration does not involve much more paper vai'k, 
since the coefficient a^ had been secured. It is recommended 

to use the following expression which is readily obtained by (2a): 


1_ 

2k 



-• b 'f k 
2k i 


(8d) 


Thus with (8a) to (8d) the coefficients of the power series (l) 
can be computed without excessive paper work, when f and Uq are 

given as functions of x. As a rule f is the given quantity, 
while Uq is the factor looked for. But befox-e proceeding to the 

solution of this problem, simpler formulas for the first approxi- 
mation of the coefficients are indicated. 


(b) First -Approximation for a 2 , aj + , b 1? and bo 


If the principal flow -direction is given by the u component 
of the velocity, as is usually the case for nozzles, and which 
always holds for the throat of a converging-diverging nozzle 
particularly, the value of the velocity u B , obtained by simple 
flow filament theory, represents a first approximation for the 


velocity at the axis. 


Compared to Uq the terms 



small by assumption and a satisfactory first approximation is 
secured for these quantivieB when the velocity u Q and itB 

derivatives in (8) are replaced, by the corresponding quantities 

Ug. 


are 


of 


Bearing in mind that f = 1 at the thpoat. the equation of 
continuity from which the flow filament speed can be deter min ed, 
reads 


( 10 ) 
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Differentiation with respect to x readily yields then the velocity 
increment 

u »f = - — f * (10a) 

s a 

H us 

At the throat the equation he comes indeterminate. The application 
of L’Hopital^s law 

f ! = Oj u f = — i— \/( it - l)ff** (lOh) 

s K + 1 * 

In this equation f could of course he emitted, since it was 
taken equal to unity, hut, in view of the repeated appearance of 
the expressions u’f and ffM, the form is retained. Quantity ff n 
can he 'interpreted as the ratio of half the nozzle height and 
curvature radius of the nozzle edge at the throat. 

Thus the insertion of the flow filament values in equation (8a) 
gives the following formula for b-j_f (+ sign and several dashes 
indicate that a first approximation is involved ) 

h^f = Ugf* + . . » (lla) 

By (2a) 


a 2 f 2 = u s ff * 1 + u B ’ff ! - u s f l2 (llh) 

\ 

The last equation could have been obtained .Just as well by 
the insertion of the flow filament values in equation (8b) . ’With 
the aid of (10) , u s ,. and u 0 'f, may be regarded as defined, equa- 
tion (ll) contains no unknown factors. After using the flow filament 
values for computing these coefficients instead of the axes values, 

It is logical to employ them also for computing bg. 

The boundary condition (9) Is written in first approximation 

v f = u 8 f T . . '• 

which by means of (lla) gives 

\ b^f 3 = u_f 1 - b n f = 0 

6 -3 8 x 


(llc) 



8 


NAOA TM No. 1215 


hence, owing to (2a) 


1 _ 

2k 



(lid) 


The last two -results are to he expected beforehand. 'When inserting 
the flow filament values in the coefficients in (l) it is not to 
he expected to- immediately obtain two terms of the development. 

Nothing is said further about (lla). The identical result 
could be obtained without any theory by an appraisal of v with 
the aid of (9). Formula (lib) is more instructive. Forming u 
with its aid, w can be calculated in first approximation. 


w = u 

~ 1 f Y 

1 + *r ( — 

.f + . . .1 

= Urt ‘ 

1 + i 

J 2 

/ u *f V 

(ff tf + B f*J 


2 

j j 

0 

2 

f 2 

\ % /_ 


(lie) 

Interesting moat of all is the last bracketed term; quantity ff* 
is always positive at the throat. ••• Assuming, subsonic- velocity in 
the converging part of the nozzle and supersonic velocity in the 
diverging part, the last Summand in the subsonic zone i3 negative, 
in the supersonic positive. Two effects can be differentiated in 
first approximation, which cause a deviation in velocity w at a 
point y from the respective velocity Uq. One of the effects 

rises as a result of the curvature of the nozzle edge, and is given 
by the term ff tT , >the other is caused by the inclination of the 

u a 'f 

nozzle edge and is given by the term f*. In the subsonic 

u 

s 

range, curvature and Inclination effect carry opposite signs; in 
the supersonic range the effects are accumulative. So if curves 
of constant velocity are plotted in the nozzle, the lines in the 
supersonic range will be curved much more than in the subsonic 
range. In fact the curves of constant velocity may even run 
parallel to the y-axls in certain circumstances. This character- 
istic behavior of curves of constant velocity or constant pressure, 
which is the same, are continuously observed again. 
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(c) Derivation of a Differential Equation 
for u Q (x) "Volume Corrections 

; 

To complete the solution the function u Q (x) must "be found. 

The character of the profiles is determined by the coefficients. 

It is clear that an equation for u Q must "be obtained with the aid 

of the continuity condition, •which is written in the form 

tb JL 

u(l - u 2 - v 2 ) K 1 dy = M (12) 

Jo 

M is a mass flow made dimensionless by the tanfc density, the 
maximum velocity and the throat cross section. 

An equation for Ug is obtainable then by developing the 
integrand at point u « u Q and v = 0 up to the terms of the 

fourth power of y followed by integrating. But a somewhat 
shorter process is preferred which in addition brings out the 
importance cf the individual terms more clearly. 

In analogy to ( 7 ) 

1 

e(u, v 2 ) = u(l - u 2 - v 2 ) 1 * 1 (7a) 

Observing that the development of u - u.^ and of v 2 starts with 

the second power of y, the development for 6 (u, v 2 ) can be 
confined to 
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^ iS J cea ^--^S r aeon that this equation can "be expressed differently 
TChen function ( 7 ) is introduced everywhere: 

e(u, t 2) = a 0 + e ua ( u - % ) 



Introducing in this expression the development of u - uq and of y2, 

vhore only terms up to including the fourth po^r of y a ^ e con- 

S ' Si perf ^ n f the int 0 S ral 111 U 2 ). the application 
of (oa), ( 8 b), and ( 8 c) gives 
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f *u. 


+ -|i-V -li- * 2 t\t 


2 * V 5 u 0 


5 u r 


1 Ij. 1 -3 

+ -— at if 1 + — b.f- 5 

24- ^ 30 3 


4 


+ 9 ( a 2 f ^ + — — a p f 2 b 2 f 2 

uuo\ 40 2 20 u Q 2 1 


0 


40u 


uuo 


- 3 


uo 


U 


+ 3 


0 


®0 


u. 


4-> 


H>i f 


This equation can be considerably simplified when certain 
limiting assumptions are made. Bearing in mind that the use of 

1 * r k 

Hi. T 

120 


still other terms in equation (l) besides a^f" 1 ’ "would also 


1 ■ 6 

involve a term of the form — - agf it is apparent that there is 

7 • 

no sense in carrying summands amounting to a mere tenth part 

of a, f 1 * or — bof3. The same holds when they are (smaller 

120 4 ^0 3 _ 


1 P 

than the hundredth of -g- agf . 

To effect a substantial simplification consistent -with accurate 
results the task is restricted to an area near the throat section. 
For this area it is assumed that 


and 


if* I < 0.10 


ff M < o .50 


(13) 


(1*0 


With the aid of (10b) and (lla) it immediately affords, the folio-wing 
estimates for k = 1.400, that is, say, for air 
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Since uq has at the throaty in general, the value of 0.40, we 
get 


6 u 


~ (V,f + u ’f) S 0.1 

ill ' X 0 / 


I -~~ o v) " °- 01 */ 


All terms in the expression in the brackets can he voided except 
the underscored term; the terms are small enough so that there is 
no danger that added up they could contribute a substantial amount. 

In the last bracketed expression the greatest term is given by 0 UUO , 
as seen in table I. This term cannot exceed the value 15. Since b-^f 

occurs in the fourth power there is no hesitation to omit the last 
term also. The result is the following equation to which on the 
assumption of (13) and (it) very substantial accuracy is attached; 



+ 


uuo 




(12a) 


Equation (12a) may be written in the form 
M 


f = 0 O + 9 uo 


+ e 


[j a 2 f2 + 35 + ? f ' ( b l f + V f )] 

e u«o [j a 2 f2 + 120 a l/ + ? f ' ( l l f + 
uo [90 " 


— a 2 f 2 a h ? k 
720 d 4 


2 /* -p - p2 


+ 20^ a 2 fC (V) - IS “2 


f 2 (b x f + u ¥) f* + . . . 


The development of 9 wa s restricted to terms of the second 
derivative at the most, which is ample for the smell velocity 
differences. In the last equation the first three terms can equally 
be regarded as development for the value 
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u 0 + du = u o + I a 2 f2 + Ti5 + 6 f r (V + u o ,f ) 


(15) 


at point u = u Q , hence ■write the equation in the form 


M = f0(u Q + d.u) + f 9 


uuo 


I 


v2 


720 




V , , » 


( 121 ) 


The brackets for 0 contain the value of the function argument. 
The second s umman d contains only the terms related to the throat. 
The principal term is the first one* as seen from (l) . 


Since 0 UUO is invariably negative, the last term represents 

a small negative quantity, which is easily estimated. At ff" = 0.20 
it corresponds to a 1:5 ratio of half the nozzle height to curvature 
radius at the nozzle edge of the throat; a first approximation 

6 

by (lib) gives a of 2 =• 0.03; ~~ is less than 15. The second 

D 

summand of (12b) is in this instance, not quite 0.001 part of the 
first. A later accurate calculation gives it at 0.5 percent of the first, 

l’ f 

The through flow volume is nothing else but j 9 dy at the throat, 

CO 

but being near the maximum of function 9, the value of the 
integral of is but very little different. Consequently, 

the through flow volume of a Laval nozzle should approach, as a 
rule, «the value from flow filament theory very closely. So, if .the 
portion of the through flow volume by which this is smaller than by 
the simple flow filament theory is designated as volume correction dM, 
equation (12b) can be written with the aid of (10) as 


0 s f + dM = 


f0 (u Q + du) + f0 uuo 



— a^f 2 "- ^ 

720 


aof^ai.f ' 



(12c) 


The principal terms here are the first term at the right- and 
the left-hand side. Since the last term is very small, its 
variation with progressing x near the throat is very small, too, 
hence a solution by equating the two principal terms. From 


e s = e ( u o + 4u ) 
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the equality of the arguments can, of course, not he deduced, since 
the reverse function of 8 is ambiguous. However, it is obvious 
that to ,a subsonic u s a subsonic value uq + du must correspond 

and the some holds for supersonic value, else the difference 
of u s and uq + du would continue to increase with increasing 

distance from x = 0, whereas the two values are certainly near 
each other. Hence by (15) 


r 2 + — 


.The relative volume correction is then 


to 


the values at the throat to be inserted at the right-hand side. It 
is of insignificant influence when the following approximation 
formula is used: 


0 uuo 9 uus t (k + l) 2 

'2 ' ' ' “ "T *»" • * • 5=2 “ 1 _ '!“••• 

®max ®m ax K - 1 

Equation (l 6) has a special characteristic. Select a place 
at the. nozzle where f f = 0. It need not be the narrowest place 
of the nozzle, and average the u component of the velocity over 
the cross section. This average may even be equated to the 
average velocity w. Although v is in no ways equal to zero 
within the limits of accuracy for f* - 0, v still is small 
enough that it plays no part for the formation of w, where it 
enters squared. Hence for f* = 0 




On considering equation (l6) it, is seen that at the throat 
or any other place with parallel wails the velocity averaged over 
the section is equal to the flow filament velocity u s . 

For the derivation of (16) the assumption had been made that 
the relation 
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1 


aM = jp e 


mo 


I 


1 

720 


a 2 f 2 a^ + 


is satisfied; tut this is complied with exactly -only a; the throat 
■where the correctness of the equation was simply demanded. 


To ascerta in the extent to ■will ch the fact that this equation 
does not hold at other points of the nozzle may affect the equa- 
tions (l6), we put 


dM - f 0 


uuo 



1 

720 



+ 


ef 


Uq + du is then not exactly eq.ua 1 to Ug and differs, say 
hy dUg, from this quantity. 


u Q + du = u B + du s 

Introduction of ef in (12c) followed by development gives 


9 b + e = 9 ( u a + 4u s ) » 0 
that is a quadratic equation for 


s + ®us^ u s + g ®uus ( du s) 
dUg with the result 


2 



But the practical calculation shows that the second summand under 
the root is small even very close to the throat with respect to 
unity. Therefore the development gives 


f* 4 0; 



]Por a nozzle with the value ff ,T = 0,20 at the throat the 
following relationship between the error dug and the value f * 
was ascertained: The bracketed expression contains a few more 

terms that were omitted in the equation of ef . Besides, the 



1 6 


MCA EM Uo. 1215 


value was computed also within a region of f * which no longer 
satisfies (13). Figure 2 indicates that the error is entirely 
insignificant. It amounts to less than 0.5 percent of the velocity. 
So, within the limits of error of the theory, equation (l6) can 
he regarded as exact solution of (12) and (12a) . 


(d) First Approximation of Uq and of the Yolume Correction, 
Range of Yelocity of the First Approximation 


As in the : foregoing a first approximation for Uq is secured 
by the Insertion of the first approximation of a 2 f 2 , a^f^, 

and b^f in equation (l6) . Quantity Uq ' is, of course, replaced 
by u^ 1 


Uq = u a 


1 



ff tT +2 



(l&) 


Thus the velocity distribution in a nozzle can be computed by (lla) 
and (l6a) when i.ts principal flow direction is given by the velocity 
component x; u a and 'f are obtained by (10), (10a), and (10b). 

Equations (l6) and (l6a) are significant by themselves. In many 
cases the pressure distribution is measured in the nozzle axis. 

To compare it with the theoretical results quantity u n can be 

computed by (l6) or (l6a) depending upon the degree of accuracy 
required and the pressure at the axis readily obtained. If 
satisfied with a first approximation, a direct relation between 
the pressure at the axis and the pressure of the flow filament 
theory is simply derived according to (l6a) . 

Equation (l6a) resembles (lie) very much. Forming, for 
example, the average velocity in first approximation, we can, 
because the expression 


laf 

2 f2 



+ 


u 'f 
% 
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must be small relative to unity, ■write "by (He) 



where, of course, Ug - w = u for f * = 0. On the other hand, it 
is seen that this equation is not valid for f * ^ 0 since the term 
■with f* in (l6a) contains the factor 2, out the qualitative 
statements regarding curvature and inclination effect remain 
applicable as before. "While in the supersonic range u Q must 

always be smaller than u s , u 0 can be > or <_Ug. 

Equation (l6a) also affords a means for the range of validity 
of the first approximation. In the coefficients snd b^f the 

quantity u v a can be replaced with more justification by u s 
as the slope f * of the nozzle and the dimensionless 
curvature ff 1 * become less. Forming with (lba) the first approxi- 
mation of u^'f gives 


u r Q f = u a 'f f- |(u s f 2 ft< * + 3 u s 'f ff»* + u B f 1 ff" 

+ 2 u_*f f* 2 + 2 u ” f 2 f‘) + . . . (IS) 

B B 

From this equation the question of substitution u Q 'f by u g f 
Is seen to be much more difficult to answer. At any rate it may 
be stated that this cSn be done with less justification at the 
throat of nozzles symmetrical about the y-axis than for Uq 

and Ug, for by (l8) 



u 'f 
u 0 




where ! f j 2 also enters in the formation of a 2 f 2 , a n d In it 

the main reason for a persistent failure of the first approximation 
Is to be found. It must be reckoned with that for ff * * small 
with respect to unity the error of the first approximation 
near f* = 0 is of the order of magnitude of ff 1 1 itself, which 
is, that the first approximation differs too much from the flow 
filament theory. Therefore, the applicability of the first 
approximation Is predicated on the assumption that in the development 
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of 1 + ff M with respect to ■ ff»* the series can he broken off 
after the first term. 

Equation (lj) affords the relative volume correction factor, 
that is, the difference of actual through flow volume from that 
computed by flow filament theory. Putting in uke equation the 
first approximation results in 

M = _ iLXi'(ff »*) 2 + . . . (17a) 

Mg 90 

Considerable inaccuracy must be presumed for. greater ff M , because 

the same statement made about the calculation of u^, where agf 
was employed, applies in a considerably .greater measure here, where 
this quantity appears squared and a^f^ appears. Subsequent 

examples will show that this formula exhibits very substantial 
errors, where the other first approximations still give very good 
service. But it is practical as estimation formula. The volume 
correction for nozzles to which the present theory is applied, is 
unimportant . 


(e) Solution of the Equation for Uq 

What is the best solution of (l6)? Since agf 2 , b^f, 

and a^f 4 are in the final analysis dependent upon u Q and its 
derivatives, an ordinary differential equation for uq is involved. 
By (8a) and (8b) 


u- a u 
0 s 


■— - - 

120 4 6 




(u ’f) 2 + u «f f * - u ~~ f — (u ’f) 
Vo J o o e Q dx v o /] 


( 1 6 b) 


a^f 4 has been left, for the rest the equation contains only 


functions of u Q , such as G n and 


0 


7 uo 


u 


derivatives of u^ and 


functions of x, such es f . and u s . The greatest terms are u Q 

and u„, while — - at.f 1 ^ represents the smallest term, when the 
0 120 4 

bracketed expression is regarded as ono term. 
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Bearing in mind that -i- at-f^ is only one-fifth of the value 

120 

which the last term in the development of u (equation (l)) can 
assume/ it is clear that the quantity uq, given "by 


Ul = u o + ife a ^ fk (19) 

differs very little from Uq. Quantity uq is naturally a much 
better solution of u Q than -what "was called "first approximation" 
of Uq. - The derivatives of Uq are therefore very closely 
reproduced by the derivitives of u^_ • and for the calculation 

O 

of a 2 f and bqf, -which themselves represent only portions 
of Uq, quantity uq is certainly sufficient. So if this 

quantity is used for computing those two coefficients equation (l6) 
can be written as follows: 



+ Uq'f f * 



(l6c) 


It is a differential equation of the second degree and second 
order for uq with very characteristics! properties, as will be 

shown. Quantity Uq can now be determined by an iteration process, 
approximate values for uq and its derivatives being entered at 
the right-hand side of the equation and a new u-^x) calculated, 

This method however does not always give a very quick result. It 
is less disturbing to make a double graphical differentiation 
of uq because of the last term. In view of the smallness of the 

particular term near sonic velocity the error introduced here is, 
in general, unimportant. Much mere disturbing is the fact that 
the error in u^'f may substantially exceed that in Uq. Naturally 

the iteration is best started with the insertion of the approxi- 
mation by (l6a) and (l8), which can be just as well regarded as 
first approximations for u 0 and UqT as for. uq and u^'f. 

Figure 3 shows the convergence of this. method for the nozzlo a = 0.20 
(described in BoctionIVc) . Its value at the throat x = 0 is 
ff lt = 0.20. 

The x-axis represents the result of the flow filamont theory, 
a subscript added for uq indicates the number of iterations; 
u r is the first approximation, obtained direct from the formulas. 
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The second approximation is not as yet quite satisfactory, hut 
even so it lies in vicinity of the throat x = 0 very close to 
the solid curve -which represents the final solution. This -was 
obtained by means of a method -which yielded several solutions in 
the subsonic range x< 0. The solution plotted here has the 
property of having a common intersection point about equal to -1.2 
with all approximations for x. Three to four steps are necessary 
to reach a result corresponding to the accuracy of the method. The 
convergence of this iteration is likely to be largely dependent 
upon the nozzle shape. In many cases where the first approximation 
is still a little uncertain, the second should give ample service. 


Another method of solving (l6c) gives quicker results. As 
mentioned in the foregoing, it is seen that solutions are lost in 
the subsonic zone with the iteration method. We shall deal with 
these cases where the flow in vicinity of the throat rises to 
supersonic speeds. After reaching sonic velocity ©^p = 0 the 

last term plays a subordinate part in the neighborhood of the 
throat. It is advisable to first establish the first approxi- 
mation of up and Up 'f which are nothing else but the first 
approximation of u Q and u Q 'f, by means of (3.6a) and (l8). 

The magnitude of the last term can then be easily estimated by a 


second differentiation, since the functions 


G- and 


U Q 


u 


are 


tabulated. In the neighborhood in which the last term is only 
about 10 percent of the right hand of the equation the equation 
is best regarded os differential equation of the first order in 
and then solved with respect to Up'f. 


U 1 


-f f + 


ui ,f “ 


° 1 j 6 ( u = - u i) + f fj- (V f ) U 1 tO 


I 


+ f»2 


G-i 


(l6e) 


The characteristics of this equation are best seen when the small 


terms f*. 


and 




are regarded as nonoxi3tant. 


If Up is chosen too small, the value for Up' will be too high. 

Hence in a progressive integration in positive x direction, that- 
is, in flow direction, a3JL integral curves in throat vicinity lead 
asymptotically to the same end curve, independent of the initial 
value. If the size of the steps is not chosen small enough a 
strongly damped oscillation of the values result instead of a one- 
sided approach to the asymptote. Therefore a solution of (l6e) is 



MCA TM No. 1215 


21 


secured by inserting close approximation values under the root 
and then defining by progressive integration in flow direction. 

"While for u^ only an initial value must be assumed, — (u-^'f) 

must be taken from the first approximation until the curve of u*pf 
is so steady that this curve can also be differentiated. Of course, 
the differential quotient of u^'f can never be computed at the 

exact place where 'f itself is looked for, but this iB unimportant 

at the smallness of the particular term. It is practical to plot 
the. departure from the stream filament velocity Ug - up and u^'f 
against x in the interpretation of (l6e) and to plot the curves 
of the stream filament theory, and of the first approximation before 
starting the calculation. In the first case the flow filament theory 
naturally gives the x-axis. 


If the initial value of up is chosen incorrectly, which is 
evidenced by considerable fluctuation of u-> ’f during many steps, 

it is better to start all over again, in order to obtain a smooth 
curve before the throat as soon as possible. Equation (l6e) is 
very appropriate even in the entire supersonic range. 

The described behavior of (l6e) thus indicates that values 
of u-j_ and u-^ l f are obtained at the throat which are practically 

independent of the chosen initial conditions. Physically it means 
that there is always a tendency toward a well defined flow attitude 
at the throat. 

In the range in which the last term of (l6e) already playB an 
essential part, that is, exceedB about 10 percent of the right-hand 
side] it is better to apply a formula obtained by integrating the 
just cited equation. 


Up*f 


rtx e f 1 

* f ? G l(V f ) 2 + f 'V f - - u l) ^ + ( U l' f ) 


(l6d) 


Having defined by (l6e) the piece' of the' curve of u^¥ and u-j_ 

near the throat, the whole curve in the subsonic. and supersonic 
range can be computed by (l6d) j Xq is tho point at which the last 
calculation is started. Approximate values of u^ 'f therefore 
substitute for x in the integral, as obtained by extrapolating 
the u^f curve. The determination of u^ f f proceeds in the 

subsonic range in the direction toward the floor and it is found 
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that the Integral curves In this direction spread out and after 
several steps lead to perceptibly different results. It results 
in deviations from the curves obtained by iteration, in both direc- 
tions, while in the supersonic range and in the transitional zone 
all the computed curves are coincident (fig. 8) . The next probably 
surprising fact of spreading of the Integral curves in the subsonic 
range has a very good sense physically. Flows, which in the sub- 
sonic zone exhibit marked differences, still give the same flow In 
the narrowest section of the nozzle, This agrees with the experi- 
mentally known fact that the flow in supersonic nozzles iB almost 
Independent of the shape of the nozzle entry. It probably alt/ays 
tends toward the state of flow for maximum through flow volume. 

Attempts to start the integration in flow direction for 
arbitrarily specified initial values up and u-jf in the subsonic 

range do not lead to the desired result. Even a slight variation 
of u-^. at greater distance from the throat for constant u, 'f 

produces . prof ound deviations before the throat. Figure 4 illustrates 
two examples (the same nozzle is in fig. 3). 

The solid: curve Is the solution for general passage through 
the sonic velocity near the throat. The choice of a slightly too 
high value of Uj_ (dashed curve) at x -• -1.5 results in a climb 
of u^'f beyond all limits. The flow cannot- be continued at all 

through the narrowest cross soction. The choice of a little too 
small a value of up (dash-dots) leads to solutions at which the 

velocity of sound Is not reached at all.- Near the throat u, r f = 0 
and is ultimately negative. 

The previously cited asymptotic movement of the values to one 
and the same end curve therefore refers only to a very narrow range 
near f l = 0. 

It is seen- that it Is now possible to define up and Uj_'f 
as function of x. The coefficients a 2 , a^, b^, and bo can . 
be computed by (8) where Up Is utilized Instead of Uq with 
sufficient accuracy. Equation (19) then yeilds Uq and equa- 
tion (l) the velocity components u and v. 
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IH. UEETVAD3I0N OF THE FIRST APEROHMATION FOE THE 
CASE OF AXIALLY SYMMETRICAL FLOW 


In general the slope and curvature of the -walls in axially 
symmetrical nozzles is not as great as for the two-dimensional 
problem, hence the first approximation .is here of greater interest. 

Bather than repeat the previous derivations, a new method is 
employed. With x denoting the nozzle axis,- y the radial 
coordinate and f the nozzle radius the equation of continuity is 
written in the form 

(pvy) « - 7 ~ (pu) (20) 

where u and v are the velocity components in x and y direc- 
tion; and., as before, w 2 = u 2 + v 2 . 

Now, in order to obtain an approximate, value for v it must 
suffice in the above equation to insert an approximate value for pu, 
since v is to be small only relative to u. Eestricted to the 
value of the stream filament, designated by p g u g , its equation 

of continuity is 


f StP s u s = M (20a) 

As before f is made nondimensional by the radius at the throat 
of the Laval nozzle, p by the chamber density, p„ end u by the 
maximum speed, so that M is a, through flow volume which is made 
nondimexiaional in exactly the same way as in (12) . 

By logarithmic differentiation of (20a) 


_1_ d fo sUs ) _ 2 _ ?us u , 

Ps*s ^ ~ f _ e S US 


(20b) 


It should be noted that (10a) and (10b) are valid only when f is 
proportional to the cross section. If the cross section is 
proportional to the quantity f 2 as in the axially symmetrical case 
the factor 2 is additive to f* in (lOu) and to ff fI in (10b), 
as equation (20b) also indicates. 
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By means of the last equation 

^ f ^ 

( P ^ r ) - ^,,”8 f + ' * * 

This is, of course, an approximate equation, because the relation- 
ship of product pu and y is not considered at all, but an 
approximate value introduced for it. The right-hand side of the 
equation merely contains functions of x, aside from y, hence 
integration with due regard to the boundary condition yields 

y = 0; v = 0; v = b-^y = u s f T ~ + . . . (21a) 


■with — = 1 on the right-hand side according to the earlier 
P s 

omissions. With observance of the freedom from rotation we 
immediately get 

U = u 0 + ~ a 2 y 2 = u Q + ~ (u s ff »* + u s ’ff ‘ - u s f* 2 ) + . . . (21b) 


The formula for u and v agrees in first approximation with the 
formulas for the velocity components in the two-dimensional case. 
Quantity uq is again computed by the continuity equation in the 
form 


M = 2 rtf yu (l - u 2 - v 2 ^) dy = 2ix 

if 0 

For d(uj v 2 ) the Bame development as before is used (174) , but 
the result after integration is slightly different, namely. 


ir 


ye 


(u, V 2 ) 


ay 


( 22 ) 
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M 


- ** + e - (t a ^ f2 + s k 6 ^ * I y + 


+ e 


mio 


h C a 2 f2 ) 2 + • • ■ 


= 0ju n +.. r- s 0 f 2 + b n 2 f 2 + ~u 


l 0 ' 4 ^ " 4u 0 U 1 


8u 0 ^*V 


\ 

! 


+ 0 


uuo 


5e + • • ‘ 


« 


The omitted terms play no part in the first approximation. In the 
last equation the "brackets for 0 contain the argument; this is 
no product* The solution is found the same as "before 


u - u fl = r/a 2 f 2 + i- b 2 f 2 + -i-u. , fb 1 f\+ . . . 
s Q k l * • u Q x 2u 0 0 x j 

[u s ff * * -2- u„ 'f f *\ + . . . 


1 

k 


2 a 



ThiB equation therefore differs in several factors from (l6a ) , the 
corresponding first approximation of the two-dimensional problem. 
The first approximation for the relative volume correction follows 
as 


aM 


K-+ 1 

96 


(ff**) 2 + . 


Its difference from (17a) is very little. 


(sfc) 


1 7. APPLICATIONS AND MODEL PR03EEMS 
(a) The Two-Dimensional Source for Compressible Plow 


The formulas of the first approximation for two-dimensional 
flow are checked against an exact solution. Chosen is the example 
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of compressible source flow, by which a flow between two flat walls 
sloped at angle 2a is given. The problem can be solved generally, 
when assuming that angle a is small. Consider the flow at point 
of the "nozzle” of cross section f . The point of intersection of 
this cross section with the x-axis is distant x Q from the origin. 

In addition r 2 = f 2 + x~ (fig. 5) • 

f 

Therefore f*= tan a = — ; f f r ' = 0. The angle a is measured 

x o 

in radians, so that the stream filament velocity at the point of f 
is given by 

fe ( u s) - “o 9 ( u o) 


This equation can be developed. 


Putting 


a - tan a - 


•i tan^a 
3 


gives 



With observance of (10a.) and omission of terms of higher order 



1 2 i 

u - i -2- f * u - u + ± u 'ff » = 0 
3 30 us 0 B 3 S 


or exactly the equation ,(l6a) for f f ' 1 = 0 

The velocity w f in the point with the coordinates x « x Q 

and y = f can be expressed in simple manner. Tirst, it must be 
equal to the velocity on the x-axis at distance r from the 
origin, developed it gives 

w f - % * “o' ( r ■ *o) 

Second, w^ can be expressed by the x component of the velocity 
in this point u^ and the cosine of the angle, 


w f = ^ 
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Development of r gives 


r 1 p 1 

— = 1 + — f end r - x n = — ff ? 

x Q 2 0 2 

'which ■written in the last two equations gives after eliminating w^ 

“r - n 0 = -|v ,a + |V"' 

The same result had "been obtained for ff tf = 0 hy (lib) and (l). 


(b) Comparison of Theory and Experiment on an 
Axially Symmetrical Nozzle 

T. E. Stanton (reference 4) measured the velocity distribution 
on the axis and along a line parallel to the axis whose distance 
from the to 11 at the point of minimum section amounted to about 
2,5 millimeters, where ff** = 0.210, As explained in the two- 
dimensional case at such a high value of ff * * no complete agree- 
ment of the first approximation with the actual flow is to be 
expected. Our theoretical values on the axis are a little too low 
and on the wall a little too high. Stanton measured at various 
chamber pressures (fig. 6 , c* is the velocity of sound at M = l) , 

In the first three casos A, B, and C the chamber (or tank) pressure 
is so low that the maximum velocity is not attained at the throat 
of the nozzle, with increasing section the velocity drops again. 

In proximity of the wall the velocity of sound in caso C is even 
exceeded in some areas. Thus the problem, here Involves a solution 
which is symmetrical to the minimum section of the nozzle. This 
symmetry is, of course, somewhat disturbed by boundary-layer 
effects. The symmetrical solutions are predicated upon same 
quantity in order to be able to determine tho theoretical curves. 

In the present instance the axial velocity at throat was taken 
as specified, in consoquonce of which thoory and tost are nearly 
in complete accord on the axis. The computed velocity in proximity 
of tho wall iB, as in case E - F, generally too low, but this marked 
deviation is at any event, attributable to the method of measurement.. 

In the caso E-F the tank (or chamber) pressure is already higji 
enough so that a general transition to supersonic speeds in proximity 
of the throat is involved. There the calculation of tho speed 
necessitates no further assumption. The theoretical values of the 
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speed at the axis are a little too low, exactly as ejected in the 
first approximation. The experimental wall velocities in the 
supersonic range are peculiarly high. The cause of velocity drop 
in the neighborhood of x = 0.15 is undoubtedly introduced hy a 
compressibility shock. Stanton’s test series D is not reproduced, 
since it cannot be interpreted . without including boundary layer 
effects , 


(c) Flow Through Hyperbolic Nozzles, Two-Dimensional Problem 

Figure 7 represents the transition from subsonic to supersonic 
velocity on three nozzles. The nozzle edge is given by the function 

f = ^1 + ax 2 

This is the equation of a hyperbola. The half nozzle height at the 
throat is taken as equal to unity. The reciprocal curvature 
radius for f* = 0 is f f * = a. The nozzle forms ropresontod by 
the above function were preferred oyer those of Gortler and Taylor, 
since they, like the nozzles used in practico, hove the property of 
decreasing side curvature with increasing distance from the throat. 
Besides on this nozzle form a groat curvature of the nozzle edge 
for small x 'values goes hand in hand with a strong slope of the 
nozzle edges for great x values. Accordingly the smaller a is 
th© greater the accuracy of the theory. 

In conformity with the entire theory the results are velocity 
profiles rather than linos of constant volocity. But it is an easy 
matter to change to linos of constant velocity. They are preferable 
for the representation. 

The three nozzles have the following values of a: 

a = 0,10; a = 0,20; a = 0.30 

The linos of constant volocity of the final solution are shown as 
solid llneB, those of the first approximation in. dashes and thoso 
of the stroam filament solution, which aro straight, of course, 
in dashes and dots. The related velocities, expressed in multiples 
of the velocity of sound c* at M = 1, aro enterod at tho inter- 
section point with the related line of constant volocity. 

Ab is seen the first approximation for all' three nozzlos still 
provides fairly practical value, and gives a fairly accurate account 
of the deviations even in the case of a = 0.30. Tho error in the 
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final solution within the- region /which, satisfies the condition (13) 
j f * | < 0.10, is about 10 percent of the difference of this solution 
from' the first approximation. The error of .’the first approximation 
at the axis is, as already appraised for f f = 0, about a times 
the difference from the stream filament -value at this point. The 
result in the subsonic zone is not definite. Figure 7 shows the 
solutions obtained for u 0 by iteration. This always gives only 
one solution. 


Figures 8(a), (b), and (?) indicate this ambiguity for the nozzle 


Uq U Q 

a = 0.20. in the subsonic zone. In figure 6(a), is plotted 

u s 

against x; figure 8(b) shows the corresponding curve for u Q 'f; 


figure 8(c) represents the lines of constant velocity for the computed 
extreme cases. Thus the solid and the dashed curves represent exact 
solutions within our accuracy. 


Figure 9(a) is the first approximation for the nozzle a = 0.20 
for the case of maximum velocity at the throat and subsequent 
- decrease. The value i denotes the Mach nurnbor at the throat 
based on the stream filament solution. Thus for i =* 1.0 it givos 
solutions which are symmetrical to the throat and asymmetrical. 

At values of i near unity the Mach number 1 is oxceodod in wall 
proximity. . Figure 9(b) represents the practically not bo important 
casos of symmetrical supersonic flow. It attains its lowest 
velocities in the loast section, Tho i values are here so 
chosen that in one case each for 1 < 1.0 and i > 1,0 tho seme 
value of 6 , that is, tho same through flow volume, prevails. Tho 
result is that tills sories of figuros represents essentially all 
possible casos from the lowest to the highest velocities. Tho 
figures are arranged so that flows with oq.ua! through flow volumes 
have correspondingly equal piacos. Tho caso for i ~ O.98 in tho 
supersonic zono is not shown. The linos of constant velocity in 
the supersonic zone are fundamentally different from those in tho 
subsonic zone. Interesting is tho fact that at vory high velocities 
the velocity chsngo is largely transverse to tho flow direction, 
and small velocity variations aro already accompanied by appreciable 
changes in density.. 


(d) Flow Through Hyperbolic Nozzles, Axially Symmetrical Flew 

The samo function chosop. for tho two -dimensional problem was 
used for tho distanco of tho nozzle edge from tho axis f . Tho 
nozzle parameter was indicated by aj>. Merely the first approxi- 
mation, whose linos of constant velocity aro shown as solid linos, 
was computed; Figure 10 represents tho throo casos 
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= 0.10; a-g = 0.20 and = 0.30 

As the velocity changes with' the square of f , hence is much 
greater than in the two-dimensional case, the lines of constant 
velocity are closer together. A velocity change in proximity of 
the axis does not amount to as much as one in wall proximity, 
since the flow volumes are substantially greaber. In .accord with 
it, it is seen that in the axially symmetrical case a departure of 
the velocity from the stream filament value at the wall is com- 
pensated by a much higher velocity variation at the axis. For 
the rest, the great velocity differences on a section are less than 
in the case of two-dimensional flow. 


(e) Correction of Volume for Hyperbolic Nozzles 


With (17) and (17a) the volume correction can bo computed in 
first approximation from the hyperbolic nozzles described in 

section c. With 

are as follows: 


© 


■j denoting tho first approximation the results 

1 A 


a 

m 

Mb 


0.10 

-0.00018 

-0.00027 

,20 

-.00053 

-.00103 

.30 

•-.00083 

-.0024 


The volume correction is extremely small. Th addition it is seen 
that formula (If) is moroly practical for estimating. At a -- 0,30 
the value obtained in first approximation, is throe -timos too high. 


(f) The Significance of tho Wall Curvaturo at Throat in tho 
Construction of Supersonic Nozzles with Parallel Jot 

The construction of supersonic nozzles with parallol flow for 
supersonic wind tunnels 'by tho Prandtl-Busomann method generally 
proceods from the assumption that in this case. with sufficient 
accuracy at the throat of the nozzlo tho Mach number over tho 
ontiro section can be put equal to unity. To analyze the error 
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introduced by this assumption two constructions of a supersonic 
nozzle with a parallel Jet with a Mach number of about equal to 2 
were effected, one on the basis of an assumed constant. -velocity 
over the section at the ‘throat, the other on the basis of a 
velocity distribution attained by . our. method. At the throat ff" 
was taken as ff" = 0.30. Minimum directional changes' of One 
radian were involved. The volume correction of the nozzle amounts 
to about 0.8 percent. The end sections obtained with the two 
constructions should therefore differ by ‘no more than 0.8 percent 
but the difference of the end sections did amount to 2.4 percent 
and for that reason this quantity can be 'regarded as measure for 
the accuracy of the construction. The greatest sectional difference 
in the two constructions amounted to 3.3 percent so that the two 
curves may be said to show no difference within the accuracy of 
construction. The justification of proceeding with the c one true— 
tion of supersonic nozzles with parallel jet from the assumption 
Ma = 1 in the narrowest section of the nozzle is therefore 
confirmed. • 


Translation by «T. Vanier 
• National Advisory Committee 
for Aeronautics 
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Figure 1.- Section of converging-diverging nozzle. 



Sketch: The function 8 (u) for « = 1.400. 
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Figure 2.- The error du g as a function of the wall slope f T . 
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Figure 7.- Lines of constant velocity of the stream filament theory, 
of the first approximation and the final solution (two-dimensional 
problem). 
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Figure 9a.- Lines of constant velocity 
(1. appro xim ation) for s ymm etrical sub 
sonic flow. 


Figure 9b.- Lines of constant velocity 
(1. approximation) for symmetrical 
supersonic flow. 
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Figure 10. J Lines of constant velocity of the stream filament theory and the 
1. approximation (axially symmetrical problem). 


